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Abstract 

We propose type IIA supergravity solutions dual to the 1/2 BPS vacua of the BMN 
matrix model. These dual solutions are analyzed using the Polchinski-Strassler method 
and have brane configurations of concentric shells of D2 branes (or NS5 branes) with 
various radii and DO charge. These branes can be viewed as polarized from N DO branes 
by a transverse R-R magnetic 6-form fiux and an NS-NS 3-form fiux. In the region far 
from branes, the solutions reduce to perturbation around the near horizon geometry of N 
DO branes, by turning on these R-R and NS-NS fiuxes, which are dual to the deformation 
of the BFSS matrix model by adding mass terms and the Myers term. The solutions with 
these additional fluxes preserve 16 supersymmetries. We also briefly discuss these fluxes in 
the possible supergravity duals of M(atrix) theories on less supersymmetric plane-waves. 
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1 Introduction 



The AdS/CFT correspondence jT] provides a remarkable method to study the physics 
of gauge theories by their dual descriptions in string theory and its low energy limit- 
supergravity. In this paper, we will utilize this correspondence to study the dual descrip- 
tions of the BMN matrix model [2] and its 1/2 BPS classical vacua by llA supergravity 
solutions. This model can be considered as a 0+1 dimensional U{N) gauge theory which 
has a discrete vacuum spectrum and can serve as a relatively simple example of a theory 
with a finite number of vacua at finite N. This model has very similar vacuum struc- 
ture to that of the Af = 1* theory whose string theory dual was constructed by the 
Polchinski-Strassler solution 

The BMN matrix model [2] was conjectured to be the DLCQ of M theory on the 
11-d maximally supersymmetric plane- wave background [Sj, [HI- The action of this model 
can be obtained^either by matrix- regularization of a supermembrane action (e.g. on 
the 11-d plane- wave background, or from the quantum mechanics of DO branes on 
the background of the 11-d plane- wave compactified to 10 dimensions jH]. The classical 
solutions and quantum spectrum of this model have been extensively studied by e.g. 
[Hj, in], ^Oj- This model may also be thought of as a deformation of the BFSS matrix 
model ^l] by adding mass terms and Myers term to the Lagrangian. Due to these 
terms, the plane-wave background removes the fiat directions of the BFSS matrix model 
and the wave-functions of the DO branes no longer spread uniformly over the space but 
are instead localized around some fuzzy spheres. The 1/2 BPS classical supersymmetric 
solution describes DO branes sitting at the origin of a 6-dimensional subspace as a result 
of mass terms in these directions and in another 3-dimensional subspace, their matrix- 
coordinates obey the SU (2) commutation relations as a result of the interplay between 
mass terms and the Myers term [2], |H]: 

[X\X^] =^|e,,,X^ z,j,/c = 1,2,3, (1) 

where /z is a mass parameter. The coordinates X* {i = 1,2,3) are N x N matrices and 

therefore their solutions are in X-dimensional representations of the Lie group SU{2). 

Since for each positive integer n there is an irreducible n-dimensional representation of 

SU{2), each vacuum solution can be labeled by a partition of the integer N into positive 

^The action of the BMN matrix model can also be obtained by the dimensional reduction of d = 4, 
TV = 4 U{N) SYM on Rx keeping certain SU{2) invariant Kaluza-Klein modes on [Tl] . 



1 



integers rij, with = N, corresponding to the direct product of these rij- dimensional 

i 

irreducible representations. So the DO branes form a collection of fuzzy spheres with radii 
proportional to rii. For large rii, these fuzzy spheres can be well- approximated as round 
spheres up to a non-commutativity correction P^. 

These vacua have similar structure to those of the A/" = 1* theory which is a 
deformation of(i = 3 + l,A/' = 4 U{N) SYM by adding mass terms for the 3 chiral 
multiplets $j (z = 1, 2, 3) to the superpotential and the A/" = 4 supersymmetry is broken 
to A/" = 1. As a result, the F-term equations for the classical supersymmetric vacua yield^ 

TTl 

= -i-j=eijk<^k, i,3,k = 1,2,3. (2) 

Since $j are N x N traceless matrices, the solutions are also in A^-dimensional represen- 
tations of SU{2) and each vacuum is also labeled by a partition of the integer into 
positive integers with Yl^i — ^^e same way as in the BMN matrix model. 

i 

Due to their similarities, many aspects of the BMN matrix model and the A/" = 1* 
theory may be studied in parallel both in the field theory context and in their dual 
supergravity descriptions. They can both be studied from the point of view of deformation 
by relevant operators around an originally undeformed theory in the U.V. region. They 
both have fuzzy sphere vacua which can be interpreted in dual supergravity solutions as 
branes of higher dimensionality polarized from branes of lower dimensionality. 

The various Polchinski-Strassler solutions \^ in the string dual of the A/" = 1* theory 
have brane configurations corresponding to D3 branes polarized into various D5 branes 
(in the weak coupling regime) or NS5 branes (in the strong coupling regime) via Myers' 
dielectric effect [15j. On the near horizon geometry of D3 branes, i.e. AdS^ x S^, 
Polchinski and Strassler H] found that one can turn on additional R-R 3-form fluxes and 
NS-NS 3-form fluxes to polarize the D3 branes into D5 or NS5 branes with world-volumes 

X S"^. They found that the radii of the 5*^ of these D5 branes are proportional to the 
D3-charge that these D5 branes carry, under certain approximations. 

In this paper, we find similar physics happened in terms of DO branes. We start from 
the near horizon geometry of A^ DO branes which is dual to the BFSS matrix model [IB], 
and find that one can turn on additional transverse magnetic R-R 6-form flux and NS-NS 
3-form flux whose Hodge duals can couple to D2 or NS5 charge and thereby to cause the 
A^ DO branes to polarize into various D2 or NS5 branes. We find that the D2 branes 
^The three $i are rescaled to make the three masses equaL 
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polarized from the DO branes in the presence of these additional fluxes have S^-wrapped 
world- volumes and the radii of their S"^ are proportional to the DO-charge they carry, under 
certain approximations. We find that the general equilibrium brane configuration could 
consist of many concentric D2 branes each with its radius proportional to its DO-charge. 

We thus propose a holographic dual description of the 1/2 BPS classical vacua of 
the BMN matrix model, using dual IIA solutions with brane configurations, in the cases 
when all Ui are large. In the appropriate regimes of parameters, there is a one-to-one 
correspondence between the supergravity solutions in the bulk, where there are concentric 
branes carrying DO-charge, and the 1/2 BPS classical vacua of the BMN matrix model 
on the boundary, which are collections of fuzzy spheres. The concentric branes are either 
D2 branes or NS5 branes, in the regimes of weak or strong effective 't Hooft couplings, 
respectively, in the matrix perturbation theory of the BMN matrix model [Hj. In the D2 
brane descriptions, each dual supergravity solution corresponds to a way of dividing up 
and distributing the total DO-charge to several D2 branes each with DO-charge rij, by an 
identical partition of in the matrix model side, in terms of fuzzy spheres. On the other 
hand, in the NS5 brane descriptions, each dual supergravity solution also corresponds 
to a way of dividing up and distributing the total DO-charge to several NS5 branes, 
but by a dual partition^ of A^ in the matrix model side. These concentric branes in the 
supergravity side are holographically mapped to the fuzzy spheres in the matrix model 
side. In the large r region of the supergravity solutions (where r is the radial variable of 
the 9 spatial dimensions), the additional R-R 6-form flux and NS-NS 3-form flux are dual 
to the deformation of the BFSS matrix model by adding mass terms and the Myers term. 

These solutions of IIA when lifted up to 11 dimensions describe supergravitons polar- 
ized into M2 or M5 branes. They are giant gravitons (e.g. |T2j, JHl, PHj) each carrying 
a fraction of the total light-cone momentum. The light-cone momentum of the A^ su- 
pergravitons p+ = ^ are shared to several M2 branes each with light-cone momentum 

_ gj^^ with radius proportional to n^, in the same way as a partition of integer A^. 
In the M5 brane description, it is also a way of sharing the total light-cone momentum to 
several M5 branes but by a dual partition of integer A^ |13j . 

The main body of the paper will focus on the details of construction of the dual 

supergravity solutions in terms of D2 brane configurations by the Polchinski-Strassler 

method, valid in the regimes of weak effective 't Hooft couplings in the matrix perturbation 

•^Here, a dual partition of N is defined via switching the rows and columns of a Young tableau, see 
e.g. P-6- 
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theory of the BMN matrix model [8^ . In the next section, we study the R-R 6-form flux 
Gq and NS-NS 3- form flux H^, in the large r region, as perturbation around the near 
horizon geometry of DO branes. In section 3, we study the equilibrium radii of the S"^ 
-wrapped D2 branes with DO-charge in the solutions with general brane configurations. 
In section 4, we solve the fluxes and as being sourced by these polarized DO branes 
with D2-charge. In section 5, we study the metric and dilaton near each shell of the 
branes as well as in the large r region. In the last section, we discuss related issues to our 
results and also possible generalizations to the supergravity duals of M(atrix) theories on 
less-supersymmetric plane-waves. 



2 The R-R 6-form flux and NS-NS 3-form flux in the 
large r region 



Since the BMN matrix model can be considered as a deformation of the BFSS matrix 
model, its dual supergravity solutions, in the large r region, can be considered as pertur- 
bation around the near horizon geometry of N DO branes, which in string frame is (e.g. 



where R'^ = QOir^gsNa"^/'^. 

The fluctuations around this background we are interested in are the R-R flux F4 

and NS-NS flux H3 which are relevant to the couplings to the D2 or NS5 branes that 

can be polarized from DO branes. When these additional fluxes in large r region can be 

considered as small fluctuations around the above background (jHI), it's easy to see that 

the perturbations of the metric, dilaton and F2 are all of second order or higher in the 

fluctuations. Therefore if we neglect quantities that are of second order or higher in the 

fluctuations, we only need to turn on these additional fluxes without giving corrections 

to the background. For convenience, we can dualize the 4-form flux F4 into a transverse^ 

6-form flux Gq via Gq = Z'^^^ * F4, where * is the Hodge dual with respect to the lOd 

"'By transverse we raean that the forms such as H3 and Gg etc. have aU components transverse to the 
DO brane world-volume. 



m, my- 



-Z'^/^dt'^ + Z^/^dlcl i = l,...,9, 

gsZ^'\ C, = g;\Z''-l)dt, Z 




(3) 
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metric in Einstein frame. After some derivation in appendix A, the equations of motion 
for the transverse fluxes and Gq turned out to possess a simple form: 



dH^ = 0, (4) 

dGQ = 0, (5) 

d[Z-\Hs~gs*9Ge)] = 0, (6) 

d[Z-\*,H^-gsG,)] = 0, (7) 



where *9 is the Hodge dual in the transverse 9-d with respect to a flat 9-d metric. These 
constraints tell us that and Gq are both closed forms and Z^^lH^ — Qs *g Gq) is 
annihilated by both d and d*g in the transverse 9 dimensions. 

The solution should break the isometry 5*0(9) to 5*0(3) x 50(6), where the 50(3) is 
the isometry of the 123 subspace and the 5*0(6) is the isometry of the other 6-d transverse 
subspace. According to this isometry, in the large r, we should look for the fluxes of the 
form 

Hs = r'^iaTs + pVs), (8) 
Ge = r"(7*9T3 + 5*9^3), (9) 

where T3 = dx^ A dx'^ A dx^ is the volume form of the 123 subspace and V3 is defined as 
V3 = (i In r A 5*2, where 5*2 = ^ejjfcX* A dx^ A dx^, k = 1, 2, 3) and m, n, a, P, 7, 6 are 
constants. 

By plugging this ansatz, the set of equations ©-((Zj) admit four solutions in two pairs 
(see appendix B). Each pair consists of one non-normalizable and one normalizable solu- 
tion as follows^: 

first pair : 
second pair : 

H3 = aTs, G, = g;'a*gT3, (12) 

16 ^8 

H, = ar~'%T,--V,), G, = g^'ar'^^-- ^.T^ + - *gV,), (13) 

^The solution with the n = Q fluxes in (|12|l . when uplifted to lid, leads to the solution of the 
"superposition" of the lid gravitational wave and lid plane- wave, as described in [B] p. 20, |25| . 



- 3^3), Ge = g;'ar-'i- *9 Tg - - *9 ^3), (10) 

-3^3), G6 = (7rW-^(*9r3-3*9V^3), (11) 
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where the a in different hnes are different. In the language of AdS/CFT correspondence, 
the pair of n = —7 and = — 9 solutions^, i.e. (HU)), (fTTj) . corresponds to the operators 
of mass deformation in the matrix model side and the VEV of it respectively [21]. As 
will be shown in section 4, in the large r region, Hj, and are the superpositions of the 
n = —7 and n = —9 solutions, while the n = —7 solution is of first order in the mass 
parameter /i of the BMN matrix model, and the n = —9 solution is of third order in the 
mass parameter fi (see appendix F). 

Since our solutions are dual to 1/2 BPS classical vacua of BMN matrix model, we 
should have 16 supersymmetries in our solutions. This is one of the differences between 
the dual solutions of fuzzy sphere vacua of BMN matrix model and those of the Ai^ = 1* 
theory. In the Polchinski-Strassler case, the supersymmetry is broken from A/" = 4 to 
Af = 1 [22], [22], while in our case the solutions after turning on and Gq still preserve 
16 supersymmetries. 

We therefore explicitly solved the Killing spinor perturbatively in first order in /x, when 
we turn on the fluctuations of the n = —7 solution (|Tn|) of and Gq, which are of first 
order in /i. The Killing spinor before perturbation is the Killing spinor e^^^ in the near 
horizon geometry of DO branes, and the Killing spinor, after turning on and Gq 
that are of first order in fx, could be written perturbatively as e = e*^°^ + e'-^-', where e*-^-* is 
the perturbation of the Killing spinor, and is of first order in fi. 

We can thereby split the gravitino equations and the dilatino equations order by order 
in /i, and the equations for the first two orders in fi are^ (in string frame [2Hj, see appendix 
C): 



(14) 
(15) 



Va„$ + ^e* /^,rli) e« = -Q/73rii+ieV4)e^°\ (16) 



®Here n is the scaling dependence of the fluxes on r in large r region. 

^The letters with a slash denote the contractions of forms with gamma matrices, for example: /Fj 

1 TP T^ab rj ^ U j^ahc 1 I? j^ahcd 

2T-f^a6i , /Hg — gj-Habci : 4 — abcd^ 
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The first pair of equations (HH), (dH), i.e. the zeroth order equations, give the unper- 
turbed e^^'' which is known in the hterature (e.g. [221, EHl): 



e 



(0) ^ ^-1/8^ 



where r] is a constant spinor with definite hehcity via projection condition (l+r-F— )?7 = r]. 
Input this result into the other pair of equations (fT6|l . (fTTjl . i.e. the first order equations, 
we find that the first order correction e'-^^ can be separated into a time-dependent part 
e[^^ and a time-independent part e'^\ and the former has the same hehcity with the 
unperturbed e*-"-* while the latter has the opposite helicity to e^^-*. 
Explicit calculations are in appendix C and the results are^ 

6 = + eS^) + e(^\ (19) 



41) = ji{^r-x' - -Lr%")zi/2ri23ro^(o)^ ^21) 

where the fiuxes in first order in /i are 

H, = ^-]1Z{-T, + ^1^3), Ge = ^-g:^JiZ{-^ *9 T3 + ^ *9 ^3), (22) 

and Z 

r' 

We see from fj22j) that the leading terms of the fiuxes and Gg in the large r 
region are proportional to the total DO-charge N of the branes and are independent of 
the specific brane configurations in the small r region. Note that the equation for the 
time-dependent part e^^^ yields a constraint that the components Z^^[Himn — gs{*9GQ)imn\ 
should be constants and therefore Z~^{H^ — Qs *g Gq) need to be a constant 3- form in the 
transverse 9-d (appendix C). This is indeed the case, since it equals to —JiT^. 

So we have checked that our solutions with the fiuxes in large r region in first order in 

/X, i.e. the n = —7 solution, are supersymmetric, preserving the 16 supersymmetries. In 

^In the expressions of the KiUing spinors l|20(l . H21(l . the various indices are: The indices I, m, n denote 
1,...,9; the indices i denotes 1,2,3 and the indices a denotes 4,. ..,9; the indices of gamma matrices with a 
bar below are the gamma matrices in tangent space. Here we consistently use another parameter ]1 in all 
the expressions instead of the mass parameter /i in BMN matrix model. They are proportional to each 
other, i.e. fl — C/x, where C is a dimensionless factor of order 1, which may be figured out by comparing 
the D2 potential (in section 3) with the matrix model action. 
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section 4, we will see that the n = —9 solution is of third order in in large r region. In 
order to check the supersymmetry when superposing the n = —7 and n = —9 in large r 
region, we need to give second and third order corrections to the metric, dialton and F2. 
One difference from the original Polchinski-Strassler solutions ^ is that we still preserve 
all the supersymmetries of the unperturbed background after adding on the additional 
fluxes. This is not surprising since the T3 in the expressions of our fluxes is maximally 
symmetric under the isometry S0{3) x 5*0(6). The expressions for the fluxes H3 and Gq, 
as well as the Killing spinor in this section, are valid in large r region where the additional 
fluxes can be considered as small fluctuations compared to the background. The and 
Gq near the brane sources will be discussed in section 4. 



3 The position of D2 branes with DO charge in the 
small r region 

In last section, we have studied the fluxes in the large r region, which are dual to the 
operators of deformation in matrix model side. In this section, we will consider the small 
r region, where there are branes which are the holographic maps of the fuzzy spheres in 
the matrix model side. The general brane configurations in our solutions are concentric 
shells of Pi D2 branes each with qi DO-charge, where i denotes the ith shell, and 'Y^PiQi = 

i 

N . We will show that the radii of these D2 branes in our solution are proportional to the 
DO-charge qi that they carry, under certain approximation. 

In a general brane configuration corresponding to the partition = YlPi^iy each brane 

i 

is in an equilibrium position under the potential it feels in the presence of all branes. 
The equilibrium radius of each brane corresponds to the location of the minimum of the 
potential that the brane feels. Before we calculate the potential that each brane feels in 
a general brane configuration, we will first solve a simpler problem. We will calculate the 
potential of a probe D2 brane with DO-charge q in the background of the near horizon 
geometry of single-center DO branes, with the n = —7 additional fluxes H3 and Gq, i.e. 
fl22|) . turned on. Then we can generalize the result for the probe brane to the branes that 
are not probes, in a general brane configuration. 

The reason we can make this generalization is that there is certain configuration- 
independence in our solutions, for example, as showed in section 2, the 3-form Z^^{H^ — 
gs*'dGQ) remains a constant form at large r, independent of what the brane configurations 
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are in the small r region. In this section, we will show, under certain approximations, i.e. 
condition (j3(J|) . the brane potential is also configuration- independent. It only depends on 
the radius of the S"^ of the brane and the DO-charge of it, and doesn't care about how the 
other branes distribute, under this approximation. 

In the probe calculation, we require that the DO-charge q it carries is much smaller 
than the background DO-charge N so that it can be treated as a probe. However, we 
can relax this condition later when we make the generalization, due to the configuration- 
independence. The brane will take the shape of an S"^ embedded in 123 subspace and 
at the origin of the other six transverse dimensions. The DBI and WZ action of the D2 
brane with q units of DO-charge in the string frame is 



Sd2 = -td2 j d^o-e ^ \l - det(G'a/3 + 27ra;'jF„^) + td2 j (C3 + T2 A Ci) 



(23) 



where I-kol T2 = 2iTa' F2 — B2. We choose the static gauge that the world-volume coordi- 
nates are the same as the space-time ones, i.e. t, 6', where the angles parameterize the 
S"^. The radius of the S"^ in 123 subspace is denoted as ri. The DO-charge q of the D2 
brane is realized as a world- volume 2-form fiux F2 = ^q sin 9d9 A d(p. The B2 is given by 
the n = —7 solution ()22|1 of the in section 2. 

We will make an approximation that is similar to Polchinski-Strassler |4| that the 
dominating terms in both DBI and WZ actions are from the contribution of F2, which 
requires the conditions (see appendix D): 

47r^Q;'^ det F2 n'^a'^q'^rl q^r\ 27ra' F^p 2na'qrf qrf 

detG± W QsN ' Bo,p JiR^ QsN ' ^ ' 

where is the pull-back metric parallel to the 5*^. Under this approximation we can 
expand the square-root in the DBI action around det F2 and then the leading terms in 
the DBI and WZ actions precisely cancel (appendix D): 



-^029^^ j rfVZ^^/ V- det G,, ■ 27raVdetF2 + td2 j 2na' F2 A Ci = 0, 



(25) 



where we choose the gauge Ci = g~^Z~^dx^ and G,, is the pull-back metric parallel to 
time. This is because that the leading terms in the DBI and WZ actions, which are both 
contributed from the F2, describe the potential between DO and DO charges and it should 
be zero due to supersymmetry. The two leading terms are both large but they cancel. 
The subleading terms of the DBI and WZ actions are respectively: 



^1 /■ ,3 ^^3/4 V-detG„detGx [ ( 2ru2r\ \ 

rD2gs d aZ ' — = dt[ — , 26 

J ina Vdeti^ J \ Qsa q J 
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ro^lC,^ Idtf^-^ISM^ (27) 



S ^ - I dt^^ irl- ^Jlaqn) . (28) 



V ^9s 

where C3 = ^gj^jldx'^ A 82- These two terms in the brackets are of order rf and Jlrf 
respectively, which can be identified with the commutator term Tr (i[X*,X-']^) and the 
Myers term — Tr (^i^eij^X^X^X''^ in the matrix model Lagrangian. There should be an- 
other term of order Ji^rj coming from the second order corrections of the metric, dilation 
and Ci, which can be identified as — Tr | (|)^ (X*)^ in the matrix model Lagrangian. The 
subleading terms of the DBI and WZ actions, plus the second order corrections from the 
metric, dilation and Ci are expected to complete the potential with a perfect square, due 
to supersymmetry condition. This is not surprising since this corresponds to the perfect- 
square term — j dt ^Tr (|-^* + ie'^^^X^X^Y in the matrix model action [8,. So the third 
term could be in principle read off from the first two terms. 

By the approximation and supersymmetry condition, the action should then be 

2rD2 / 2 ^- ■ 
gsa q \ 3 

So we see that the brane potential depends on the radius of the S"^, i.e. ri, and its DO- 
charge q, and is independent of the warp-factor Z under the approximation ()24|) . The 
subleading term in DBI action is Z-independent since both Z factors from numerator 
and denominator exactly cancel (appendix D). The subleading term in WZ action is 
also Z-independent because C3 is Z- independent. And the third term should also be 
Z-independent due to supersymmetry. Therefore the potential of the brane only cares 
about its DO-charge and is independent of the brane configuration of the solution under 
this approximation. The potential of a non-probe D2 brane with DO-charge in a general 
brane configuration is thus still of the form (|28|) . There is a non-trivial equilibrium radius 
at 

ro ^ '^Jlaq, (29) 

which is proportional to the DO-charge q. This is similar to the matrix model description 
of the fuzzy spheres that each of them has a radius proportional to the size q of its block 
matrix-coordinates j2], [Hj, up to a non-commutativity correction [T2] . 

After knowing the equilibrium radius tq, we get the consistent condition for our ap- 
proximation from that 

^ > 1. (30) 
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This is a kind of scaling bound between q and QsN ^ where q is the DO-charge of the D2 
brane and is the background DO-charge^. The approximation (j3(J|) actually guarantees 
that the contributions to the DBI and WZ actions are both from the DO-charge but they 
cancel and then the remaining terms form a perfect square and is independent of the 
brane configuration and the warp-factor. 

So in our general brane configurations, suppose there are concentric shells of pi D2 
branes each with DO-charge qi (the label i denotes the ith shell), the total potential of all 
the branes can be considered as the sum of individual potentials and we have: 

where r[*^ denote the radii of the 5*^ of the branes on the zth shell and their equilibrium 
radii r^^ are therefore 

rj^ ~ ^mqu ^Pi<li = N. (32) 

i 

Note that there could be coincident pi D2 branes on the same shell if they have the same 
amount of DO-charge g^, which corresponds to some copies of the irreducible representa- 
tions of the SU{2) of the same dimension, in the matrix model side. 



4 The additional fluxes in the presence of polarized 
sources 

In section 2, we focused on the large r region and found out the additional fluxes and 

Gq, as perturbations around the near horizon geometry of DO branes. In section 3, we 

have studied the situations in the small r region and figured out the radius of each D2 

brane, in general configurations. The leading terms of these fluxes and in the large 

r region depend on the total DO-charge of all the branes and are independent of the brane 

configurations, while in the small r region these fluxes are dependent on specific brane 

configurations. The expressions for the fluxes in section 2 are not valid near the brane 

^This is very similar to the condition for the approximation in the original Polchinski-Strassler solution, 
i.e. ^ 1 m, where n was the D3-charge of the D5 brane and TV was the background D3-charge. We 
have different powers of n (or q) because the powers of the r in the warp-factor Z is different. In their 
case Z = while in our case Z — The warp-factor dilutes the background charge so the power 
dependence of n (or q) for different Z is different. 
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sources. In this section we will therefore study the behavior of H3 and Gq in the presence 
of these sources. 

For simplicity of the calculation, we study the special case when there is only a single 
shell of D2 branes with total DO-charge N. Suppose we have p coincident D2 branes each 
with DO-charge q = N/p, so the radius of the shell is rg ~ ^Jia {N/p). We study the case 
that the total D2-charge p is small, so the background metric can be approximated by the 
near horizon geometry of multi-center DO branes distributed on the S"^ with radius tq. 

The equation for *Fi will have a source term due to the D2-charge. We are interested 
in the and Gq on the background of the near horizon geometry of a shell of multi-center 
DO branes with small D2-charge turned on. We are now doing perturbation in terms of 
small parameter p. The warp-factor Z = ^ is replaced by the multi-center warp factor 
Zi in solution (jH)): 



which is the superposition of harmonic functions, corresponding to DO-charge uniformed 
distributed on an S*^ with radius tq in the 123 subspace and centered at the origin of the 
other 6-d transverse subspace. Here ri is the radius of 123 subspace and r2 is the radius 
of the other 6-d transverse subspace. Its easy to see that Zi reduces to Z = ^ at large r. 

The IIA SUGRA equations for and Gq on the background of the near horizon 
geometry of a single shell of multi-center DO branes with radius Tq and small D2-charge 
p should be (see appendix E): 

dHs = 0, (34) 
c^Ge = J7, (35) 



lOriro 



d[Z^\H^-g,*,G,)] = 0, (36) 
d[Z^\*^H,-g,G,)] = 0, (37) 

where there is a source term J7 due to D2-charge: 

J7 = 2^TD2gfp5{r^ - ro)6\r2)dn A dr2 A U5, (38) 

where 5^(r2) = S(x4^)S{x5)S{xq)S{x7)S{xs)S{xq) and cus = rf ■ dvo^S*^), where dvol(S'^) 
denotes the volume-form of an with unit radius embedded in 456789 subspace and 
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centered at origin. The source term appears in the Bianchi identity but not in the equa- 
tions for Z^^i^Hj, — Qs *9 Gq) and its 9-d Hodge dual^°. 

The equations for ^(ifa — gs *9 Gq) remain the same as in (jH)), Q, just with Z 
replaced by Zi. Furthermore, since when r goes to infinity, this form is a constant form, 
we infer that the harmonic form Z^^^H^^ — Qs *9 Gg) = ~/^^3- Although H^,Gq,Zi all 
depend on the brane configurations, the combination Z^^{H^ — Qs *qGq) is independent of 
the brane configurations. This also results in that the potential C3 is independent of the 
warp-factor Zi and the brane configurations. The difference between equations (Pj)-(|7|) 
and (jMj) - (jH7j) . besides that Z is replaced by a multi-center warp-factor Zi, is that there 
is a source term for Gg since we have introduced D2 sources on this shell. 

We can split both i^a and Gq into two pieces respectively: 

Ge = G^P + Gf\ (39) 
= HP+g,HGi''\ (40) 

where -^3^^ G^^ still satisfy the whole four equations (jH)-© with Z replaced by the multi- 

r(I) nil) 



center Zi. In large r region when expanded around tq, the leading terms of will 
reduce to the n = —7 solution in section 2. h'^\g^P are the contribution to the fluxes 
as if there were no D2 source. 

The influence of D2 source on the fluxes are mainly on Gg^^"* , whose equations are now: 

dGf^ = J7, d *9 Gi^^^ = 0. (41) 

The contribution of Gg^^'' is dominant over Gg^^ very close to the shell since it has the 
delta function as source. Since *gGQ^^ is closed, it can be written as 

>K9G^^^^ = (rf ^^iMn + rfd2hdr2) A cjs, (42) 

where uj2 = r\- dvol(S'^), and dvo^S*^) denotes the volume-form of an S*^ with unit radius 
embedded in 123 subspace and centered at origin, and di = 82 = h is a function 
of ri,r2. The solution (see appendix F for more detail) can be expressed through the 
function defined asY = r^'^dih, and we have 



y ^ j^[ClfrD29s>ig 1-1 



1 1 



(43) 



_[{n + ro)2 + r2]5/2 - roY + r^Y^/\ 

Since in our case the forms and F4 both have overaU factors of the volume- form of the S*^ in 123 
subspace due to the isometry S0{'3) x S0{6), the terms F4 A H3 and F4 A F4 are zero and drop off on 
the right sides of (|^ . 
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where the notations of the coefficients in Y are in appendix F. 

The leading terms of G^^^ and H^^^^ = *g G^^^ expanded in terms of rg, in the 
region where ri ^ rg, is precisely the n = —9 solution in section 2, and it is of third order 
in To and thereby of third order in fi in large r region (see appendix F). The contribution 
of Gg^^"* is dominant over G^^ very close to the shell of the brane since it has the delta 
function as source, while in large r the situation is reversed and G^^^ becomes dominant 
over G[I^^ instead. We see the consistency in the calculation of the fluxes in the presence 
of the sources since the solutions of H^^\ G'^P and h'^^\ Gq^^ in this section are the full 
solutions which just reduce to the n = —7 and n = —9 solutions in large r region in 
section 2 respectively. 



5 Metric and dilaton in large r region and near each 
shell 

In this section, we come to the discussion on the metric and dilaton. The situation is 
similar to Polchinski-Strassler |lj that in most regions away from the shells, the DO-charge 
dominates, and in the regions very close to the shells, the D2-charge dominates instead. 
This switch of the role of the dominance is also reflected in the change of the dominance 
between H^^\ G^^ and h'^^\ G^^^ as discussed in last section. We will discuss the metric 
and dilaton in two limiting regions in this section. One is in the large r region and the 
other is very near each shell. 

The general brane configuration in our solutions are concentric branes with various 
DO-charge and radii. In the large r region, since the DO-charge dominates, the metric, 
dilaton and F2 are very close to the near horizon geometry of multi-center DO branes with 
warp-factor Zi corresponding to the distributions of these concentric shells of DO branes. 
For the general configuration of several concentric shells of S^-wrapped branes with the 
ith shell having coincident D2 branes each with qi units of DO-charge (A^ = YlPili 

i 

and the DO-charge qi are all large and distribute uniformly on the spheres), the warp 
factor Zi in the solution of the near horizon geometry of multi-center DO branes, as the 
superposition of harmonic functions, should be 



2. = 5: 



Rl 



lOriTQ 
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where Rj = Q^n^ g s{piqi)a'^ ^"^ and ^ ^Jlaqi. 

Now we will look at the metric and dilaton near each shell of branes, say the ith 
shell. The total DO-charge of this shell is Ni = piQi and radius of this shell is Tq^ ~ 
^fla'qi. Very close to each shell, the DO-charge no longer have dominant influence over 
D2-charge since the metric parallel to the shell expand and DO-charge are diluted. We 
can approximate the metric and dilaton near the S'^, e.g. without loss of generality, near 
the point {xi,X2,X3) = (0,0,rQ^), by the solution of Pi fiat D2 branes with B2 potential 
on its spatial world- volume [201 , [21] (in string frame): 

^'5/2^5/2 ^ 1 _ _ \ Gir'^gsa^^^Pi 

DTT Qs^i Pi ^ i i / Ot ' Uj^ 



^'15/2^15/2 \l + aluly 

where the label i denotes the ith shell. Xii,Xi2 parameterize the spatial part of the D2 
branes, Ui is the energy direction away from the ith shell of branes in the transverse 
direction and constant that will be worked out in (jl^ . 

Since there is B2 field on the D2 branes in solution ()45|). the D2 branes can couple 
to Ci and there is DO-charge on the D2 branes. Suppose we are looking at regions only 
near the ith shell but not the other shells at the same time. When away from this shell 
of branes such that afuf ^ 1, i.e. i'^'^[^5 ~ 1? the above metric and dilaton ()45|1 match 
exactly with the near horizon geometry of multi-center DO branes, near this shell of branes 
{xr. Xi,X2,X3 and Xa: x^, ...,Xg): 

^^'0 Pi 1.2 , 



ds - 1^ ' dt + '^.^ 5/2 ^ / +dxj 

n2p21/2 

= 9lzr = :.3 (46) 



where rQ-* ^ ^/ia'gj, R] = GOtt^ gs{piqi)a"^ ^"^ . We used the multi-center warp-factor Zi in 
approximated near the ith shell of the brane source, and pi is the distance away from 
the ith shell: 

Zi - . (47) 

10 (r?) p! 



Pi 



[in - r^o^r + rl]'/'. (4J 
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We can define a cross-over distance pt^ as tlie distance away from tlie brane wliere 
ttiUi = 1, wliicli could cliaracterize tlie regions of influence of the D2-charge. In order 
for the match to be valid, we need when away from the shell the two approximations 
rr^Sjs ^ 1, (i.e. Pi > p'i^) and Zi ^ . ' (i-^- Pi ^ ^o^) ^"^^ both vahd. This 

10 (r^^'j pf 

requires the match to happen in the region pc^ <^ Pi <^ rQ*\ so for our approximation to 
be valid we need the parameters satisfy pc^ <^ 

Under this approximation, the dilaton and metric match exactly and we thereby found 
the relation between the parameters and variables by comparing the two solutions (j45|) . 

(gni): 

u. - Pi a -i ^ ^ ^'^ p« - a ( y^I^lll^^" ^'^ 



10 [rf^ a'^/^ -Qii^g, 



(49) 



Pi 



So the metric and dilaton near the ith shell can be written as 



r,2 r21/2 

= , —TV, (50) 



103/2 (r«)%f (pf+(p« 
and as discussed above it is valid when 



r^) / „5 \l/5 
'0 / % 



p. 



» 1, (51) 



which is just 3> 1, where qi is the DO-charge of each D2 brane on the ith shell and Ni 
is the total DO-charge of the ith shell. This is the same scaling bound condition as (jHUj) 
in section 3. 

For general brane configurations of several shells, the metric and dilaton in complete 
regions are very difficult to solve explicitly, but it's clear that in large r region they 
approach the near horizon geometry of multi-center DO branes with warp-factor Zi in 
()44|) . and near each shell they are approximated as the solutions in (jHDI)- In special cases 
when there is only one single shell of branes, the metric and dilaton may be expressed 
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approximately valid in all regions. Suppose there are p D2 branes each with DO-charge 
q on this single shell, pq = N is the total DO charge and the radius is tq ^ ^fla'q. The 
solution near the brane can be obtained from (jKn|) by identifying pi = p,qi = q- And then 
we can generalize the solution near the shell to all regions by replacing the warp-factor 
Zi ^ iQ^-2p5 , approximated near the shell as in ()47p. with the warp-factor Zi in all regions 
as in The validity of this approximation is again from (jFIjl . i.e. ^ 1. 



6 Related issues and generalizations to other plane- 
wave M(atrix) theories 

So far we have constructed the supergravity solutions dual to the 1/2 BPS concentric 
fuzzy sphere vacua of the BMN matrix model using the method of the Polchinski-Strassler 
solution, which is the string dual of the Af = 1* theory. In this last section, we discuss 
some related issues or remaining issues to our discussions above. 

Each 1/2 BPS vacuum of BMN matrix model can be represented by a Young tableau 
, P and it can be interpreted as either concentric D2 branes or concentric NS5 branes, 
in different regimes of parameters |Hj, JB]- The configurations in terms of concentric shells 
of D2 branes have their validities as dual descriptions when the effective 't Hooft coupling 
in the matrix perturbation theory of the BMN matrix model when expanding around 
each fuzzy sphere is small^^ jS], JSl- When the effective 't Hooft couplings are small, 
the interactions are smaller than the harmonic oscillator energies expanded around these 
fuzzy spheres, and the BMN matrix model can be studied perturbatively around these 
fuzzy spheres jH] . The concentric D2 brane configurations are therefore good descriptions 
of the BMN vacua when these parameters are small. For fixed partition of A^, we can 
always tune fi and R to satisfy these conditions. For fixed fi and R, it is relatively safer 
to expand around a fuzzy sphere when the numbers of the coincident spheres are smaller 
and the matrix size of the sphere is larger jH]. 

We have analyzed the descriptions of the vacua in terms of concentric D2 branes in 
the regime of weak effective 't Hooft couplings and we haven't studied the situation of the 
concentric NS5 branes in detail, which are expected to be valid in the regime of strong 

^^The effective 't Hooft coupling is pi ^ \--^ ^ 8 , where p^_^'' ~ ^ is the Hght-cone momentum of the 
M2 brane on the ith shell and pi is the number of coincident M2 branes on the ith shell. 
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effective 't Hooft couplings. The NS5 branes polarized from DO branes should also have 
equilibrium radii with these additional fluxes turned on since it can couple to the dual of 
the NS-NS 3-form flux and the 3-form potential via world-volume 3-form flux. 

In fact, there could be smooth solutions that are dual to these vacua. The smooth 
solutions and the solutions with brane configurations studied in this paper may be related 
by geometric transitions (e.g. [S2]), where branes and fluxes get replaced with each other. 

There are less supersymmetric vacua and time-dependent vacua e.g. |H], [E], jBl], 
jHZI that we have not discussed. For example, it would be interesting to understand 
such as the 1/4 BPS rotating elliptic fuzzy spheres described by jM], and the 1/4 BPS 
rotating fuzzy spheres by [3^] . There are also instanton solutions which are similar to 
the domain wall solutions in the A/" = 1* theory jHE]- It would be good to understand the 
dual descriptions of the vacua in the model that is less supersymmetric and/or non-static. 

The Polchinski-Strassler type solution has been widely generalized and applied to many 
other situations in terms of other branes (e.g. [10], [H], One can conjecture 

that the Polchinski-Strassler type solution is universal for any Dp branes which can 
be polarized to Dp-(-2 or NS5 branes in the presence of the additional R-R and NS-NS 
fluxes on the background of the near horizon geometry of Dp branes and the resulting 
solutions are dual to the mass-deformed world-volume field theory of Dp branes. In 
each such solution, we have a pair of R-R and NS-NS fluxes and this is mainly because 
we have two channels of polarizations. The R-R flux is more responsible for polarization 
to Dp+2 branes, while the NS-NS flux is more responsible for polarization to NS5 branes. 

The construction of dual supergravity descriptions to the BMN matrix model may 
be generalized to those of the M(atrix) theories on less supersymmetric plane- wave back- 
grounds |13j, IS]. For a general plane- wave matrix theory with lagrangian, e.g. 

-E/^' i^f + Im.kX^X^X' - \tVm, (52) 

i 

where /r=^TijkY^'', — X] i^l^yfc)^, there might exist supergravity duals which are 

i i,j,k 

similar to the case of the BMN matrix model. In the large r region, the dual supergravity 
solutions can also be considered as perturbations around the near horizon geometry of 
DO branes by the additional fluxes and Gq, which should also satisfy equations 
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((7j), and the form Z^^{H3, — gs *9 Ge)? where Z = should thereby also be annihilated 
by both d and d*g in the transverse 9-d. When r goes to infinity it should approach a 
constant form, so the additional fluxes should satisfy the relation in large r region as: 

Z-\H3~g,*9Ge) (53) 

where T3 = ^Tijkdx^ A dx^ A dx'', and Tijk are the coefficients of the Myers term in the 
corresponding M(atrix) theory on the general plane- wave background. Thereby one can 
conjecture that the perturbation by these mass terms and Myers terms are dual to turning 
on these additional fiuxes and Gq causing DO branes to polarize into some non-spherical 
branes. The brane configuration would be more difficult to describe than the case of the 
BMN matrix model. For example, in the M(atrix) theory on the background of T-dual of 
the IIB pp-wave lifted to lid, the BPS vacua correspond to M2 branes polarized into M5 
brane, where the M2 branes distributed on a fuzzy ellipsoid ji^. So generally, the brane 
in the small r region would take the shape that corresponds to the classical static solution 
in the corresponding M(atrix) theory and it should also equivalently be the shape of a 
probe brane in the presence of external fiuxes and Gq. 
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A IIA equations for fluxes 

The IIA equations of motion in this appendix are in Einstein frame and in the convention 
of The bosonic equations of motion are: 

d * = -^Qse-^H, A *H, + ^gl/^e''^/^F2 A *F2 + ^g'/'e^/'F, A *F„ (54) 

die'^/UF,) = g.e^^'HsA^h (55) 
d(e*/2*F4) = -gy'F.AH,, (56) 
\gsF,AF, = die'^^H. + gl/'e^/'CA^F,), (57) 
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where F2 = dCi, F4 = dC^, H-^ = dB2, F4 = F4 — Ci A H^, and the Bianchi identities are: 
dF2 = 0, dF^ = 0, dH^ = 0. The relation between the Einstein frame metric and string 

frame metric is (G^.) Epstein = 9^^^^-"^^ (G'M-)string- 

The unperturbed background is the near horizon geometry of N DO branes. In Einstein 
frame it is 

ds' = -Z-'/^'dt^ + Z'/^'dl^l z = l,...,9, 

e" = g,Z'/\ C, = g;\Z-^-l)dt, Z = ^. (58) 

The perturbed fluxes and F4 are small fluctuations in large r region. The per- 
turbed terms in the right sides of equation (|^. (|^. if non-zero, are at least of second 
order fluctuations. Therefore if neglecting the terms of second or higher orders in the 
fluctuations, the equations for the dialton and F2 still have the form: 

And therefore the other fluxes ,^4 obey: 

d{Z^/^*F4) = 0, dHs = 0. (60) 

dFi = 0, d{g;^Z~'^/^ *H3 + gsZ^/^Ci A ^^4) = 0. (61) 

From (j60|l we can define Gq = Z^^^ * F4, so that dGe = 0. We can write (j6T| in 
terms of the 9-d Hodge dual *9 in the transverse 9-d with respect to the flat metric: 
*if3 = dx'^ A (— *9 H3), F4 = dx^ A {—Z~^ *9 Gq), then from we precisely arrive 
at equation (jEI), in section 2. 

B Linearized solutions of fluxes 

In this appendix we discuss linearized solutions of additional fluxes in terms of tensor 
harmonics. Without loss of generality, we can define a T3 analogous to Polchinski-Strassler 
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|1] as follows: 



T3 = dri A = dx^ A dx'^ A dx^, (62) 

^2 = ricu2 = A cix-^' A dx\ dS2 = STg, (63) 

2 

V3 = d\nrAS2 = -^dr^ AU2 + -^dr2 A U2 (64) 

= -^dx^ A dx'^ A dx^ + — — —Eijkdx" A dx^ A dx^ , (65) 



(^,j,A; = l,2,3,a = 4,...,9) 
where the definition of a;2, uj^ are in section 4. Then 



^gTg = dr2AU5, 

9 

♦gVs = ^cir2 A 0)5 ^-c/ri A u^. (67) 

If we search for the solutions of the form in large r region according to the isometry 
50(3) X S0(6): 

= r"'{aTs + PV^), G, = g;'r^{j hTs + S Vs), (68) 

where m, n, a, /3, 7, 6 are constants. The factor in the ansatz for Gq is introduced for 
convenience in this appendix. The four equations (jH), (0), (0), © give the constraints: 

TTl 

dH3 = 0, (3 = —a. (69) 



71 

dGe = 0, 6 = —7. (70) 

n + 6 

r)2 _|_ IRr) _|_ 42 

d[Z-\Hs-9s*9G,)] = 0, lar^ 7^" = 0- (71) 

n + 6 

rf[Z-i(*9/73-^7sG6)] = 0, ^ ^^^^±^^^^ar™ - ^7r" = 0. (72) 

3 + 6 

Since in the perturbation, 0,7 cannot be both zero, we need m,n to be equal, and 
then we get 

n(n + 16)(n + 7)(r2 + 9) =0. (73) 

This leads to 4 solutions that are two pairs of non-normalizable and normalizable 
solutions in (Uni), (HH), (ini), (HSI) in section 2. 
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C Killing spinors 



The IIA susy transformation rules we used are in string frame (e.g. 

4 o o o 

= [(V9„<|. + ^e* A2rii) + (^//3r^+^eV4)]e, (75) 

where letters with a slash denote the contractions of forms with gamma matrices: /F2 = 
ft^ = j^Habc^"'"', = ^Fabcd^^'"'^ and similarly for other forms. 
The unperturbed Killing spinors e*^°^ in the absence of H^, Gq satisfy: 

[-r'"9„<i> + -e* /p2r^]e^°^ = -z~^/^r-diZ{i + r^rii)e(°) = o, (76) 

2 8 8 

[do + ^uJoabT'"'' + le" /F2rorii]e(°) = [do + Iz-'/'V^mZil + r^r^)]e'^'^ = 0, (77) 
4 o o 

[d, + ^ujiabT^' + le" /^2r,rii]e(°) = [{d, + lz-'d,z) - lz-%z{i + t^t'-^) 



+lz-'r-TidjZ{i + r^rii)]e(°) = o. (j ^ i), {7i 

8 



where the indices i,j denotes 1,...,9 and the indices of gamma matrices with a bar below 
are the gamma matrices in tangent space. So we have the already familiar result (e.g. 

m, my- 

,(0) ^ z-'/8^, (79) 

where r/ is a constant spinor satisfying (l+r-F— )// = 0. The Killing spinors in the presence 
of small fluctuations of H^, Gq can be written as 

e = e(°)+e«, (80) 

where e*-^-* is the perturbation around e*-'^-'. When we turn on the H^,Gq that are of first 
order in /x, e*^^-* is of order jj,. The variations of gravitino and dilatino of the first order in 
/i give: 

^z-^/^r^9,z(i + r^^rii)e« = -[\ fl,v^^+\e'JM'\ m 

8 4 8 

[^0 + \z~-'/'T'-T%Z{l + r^irii)]e(i) = -[li7o,,rtii + ieV^roje^, (82) 

o o o 

[(9, + ^Z~^d,Z) - ^Z-'dail + T^T^) + ^Z-'r-TldjZil + T^T^)]e^'^ 
= -[lif,,,r'^^rii+ ie^T^^r.JeW. (j ^ ^.) (83) 
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Now we can first consider the dilatino variation (jHT|) involving e^^-*. Since e^^-* is not 
time-dependent while e*-^-* is time- dependent, the time-dependent part of e^^^^ should be 
annihilated by (1 + F-F— ). So we can decompose e^^^ into two parts: e*^-"^) = e? , 

where e^^^ is time-dependent and e? is not time-dependent. 
We can then spht dilatino equation (|HTj) into two equations: 

(1 + r^^r^)e['^ = 0, (84) 

mZ{l + r^-r'-^)e? = ^3rii+leVje(°). (85) 

o 4 8 

Substituting (jHH), (jHSl) into the time-component of the gravitino variation, equation 
(I82|) . and the spatial components of gravitino variation, equation (|83|) . we have: 

doe? = ^Z-i/^F^[/f3rii- eVjeW, (9, + lz'%Z)e? = 0. (86) 

Since the right side of the first equation in ()8fij) is time- independent, we solve that e^^'* 
is linear in time: 

e? = ^ ■ l^Z-'[H,,k - gs{*9G,).,k]rm% (87) 

and e^^'' has the same helicity to e^^-*. The spatial-independence of Z^^^^e? from the 
second equation in (|HT)|l and (|7n|) imply: 

Z'^[Hijk - gs{*9G6)ijk\ = const. (88) 

The discussion so far doesn't require 5*0(3) x 5*0(6) symmetry but only that the fiuxes 
H3 and Ge be small. There are stronger constraints than merely that Z^^lH^ — Qs *g Gq] 
would be a constant from the spatial part of the gravitino variation involving e? from 



[{d, + -Z-^diZ) - ^Z-'diZ{l + F^Fii) + ^Z-^r-TldjZ{l + T^T^)]e? 
= [-lH,a,T^'' + lz'/' jG,T^]rh^'\ (89) 

o o 

where we define G3 = Qs *9 Gq. If contracting both sides of equation (j89|l with F-, and 
input equation ()85|) . and then acting on both sides the projection (1 — F-F— ), the right 
side becomes zero and we have 

(1 - T^T^)e? = 0. (90) 
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So e^2^ has the opposite chirahty with respect to e*^°) and e^^"*. 

From now on, we will use the indices l,m,n to denote 1,2,3,4,. ..,9, the indices i,j,k 
to denote 1,2,3, and the indices a,b,c to denote 4,.. .,9, for convenience. Applying the 
projection condition (IHUjl on e^2^ and then equation (jHHj) becomes: 



(1) 



3 



(91) 



Substituting (jHH) into (|H^ and using the projection condition we have 9 individ- 
ual spatial equations: 



1^; 1 



1 



(92) 



1 

12' g-™-- ■ ^ 

Now if the we look at the fluxes in the ansatz (jHJ, Q and combine the first supersym- 
metry condition from (j88j) . we have 

Hs = R'r-\-aT^ - (3V^), Gs = R'r-\^T^ - (3V:,), (93) 

where a,/3,7 are constants. 

By dimensional analysis from equation ()91|) . (j^J^ . one finds that €2^'' after extracted 
out the Z^/^ factor should be linear in so we can try the ansatz: 

e« = (/i^rv + ;U4r%'^)z3/8ri^r%7, (94) 

where /zi, /i4 are constants. Then comparing the left and right side of equation (jHH), we 
get the relation 

/^i = J7^ T7' = 77^ • (95) 
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Comparing the left and right side of equation ()92|) . we get another relation 

a + 7 a + 7 
/^i = — = (96) 

This shows /xi : /i4 = 2 : —1, and combine with and we have a : /3 : 7 = 3 : —7 : 
— 1, which is just our n = —7 solution in section 2. 



D Approximation of D2 potential 

In this appendix, we write some details in the approximation of D2 potential in section 3. 
The DBI and WZ action of the a D2 brane with q units of DO charge in the string frame 
is in dSni) 

Sd2 = -rD2 j rfVe^*-^- det(G'„^ + 27raT„/3) + td2 j {Cz + 27ra'j^2 A Ci), (97) 
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where 27ra T2 = ^ira F2 — B2. We choose the gauge that the world- volume coordinates are 
the same as the space-time ones, i.e. t, 6, ip, where the angles parameterize the 5*^. G,,, G± 
are the pull-back metrics parallel to the time and the spherical directions respectively, so 
we have: 

det G„ = -Z-^/\ det = Zr\ sin^ ^, (98) 

where ri is the radius in 123 subspace. 

The DO-charge of the D2 brane is g, so the world-volume 2-form fluxes F2 is: 

F2 = sin Qde A dip, j F2 = 2nq, (99) 
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so we have Fg^ = \q sin 6, det F2 = iF^pF'^'^ det G± = 4 sin^ 6, F^/jF"^ = ^fj^- 

Suppose the dominating terms in both DBI and WZ actions are from the contribution 
of F2, which requires the conditions (j24j) then we can expand the square-root in the DBI 
action as 



J- det(Ga/3 + ~ \/- det G„ ( 2™' ^detT^ + -—^^^^L=\ 

V \ 47ra V det F2 J 



(100) 



The leading term in the DBI part is 

-rD2g~^ j rfVZ-^/V-detG„-27raVdetF2 = - j dtdOdip (^ixa TD2g^^ Z'^ ■ ^gsin^^ . 

(101) 

The leading term in the WZ part is 

td2 J 2na'F2 ^ = j dtdOd^ (^na r^g^^ Z'^ • ^gsin^^ , (102) 

where we choose the gauge choice Ci = gJ^Z^^dx^ in section 3. The two leading terms 
from the DBI part and WZ part precisely cancel. 

The subleading terms of the DBI and WZ parts read respectively: 

^ f r. ,/4\/-detG„det G± f .r^sin^^ f , / 2TD2rf\ 

(103) 



where the Z factor cancel exactly in ()103|) and 

rD2 1^3 = J dtdOdip (^^Tn2g7'jlrlsme^ = J dt (^^^^^ ) , (104) 
where C3 = ^g~^]ldx^ A 5*2 since F4 = F4 + Ci A iJa = —gj^jldx^ A T3. 
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E Putting source 



2 



The IIA bosonic action without external source is 

Sua = ^[J d^^'x^/^R - ^ j (d<^ A *d$ + g^e-^Hs A *H, + (^^'e^^/^Fs A *F, 

2e*/2^^ A *F4 + glB2 A F4 A F4)]. (105) 

When we add D2 branes with DO-charge distributed in a single shell as in section 4, 
we actually introduced the source terms in the action. In analogy to electricity and 
magnetism, now in the action there should appear a source term (some related discussion 
on adding source term is in jUj): 

Ssource dUCs " ^2 A Ci) A J7, (106) 

where Jy = 2K'^TD2gJ^pS{ri—ro)6^{r2)driAdr2Auj5 and5^(r2) = 6{x4)6{x5)6{xQ)6{x7)6{xs)S{xg) 
The coefficient in Jy are read off from comparing the WZ action of the D2 branes with 
DO-charge involving the couplings to C3, B2. So now the total action is 

•S* = Sua + Ssource- (10"^) 

Since now C3 couples to Jj, and B2 couples to A F4 + Ci A J7, the equations for 
(i(e*/^ * F4) and for d{e''^ * H3) when adding source should be modified to 

rf(e*/2*F4) = -gl/'F,AH, + gl/'Jj, (108) 
d{e-''*H^ + gy^e''/^CiA*F^) = ^gsF^ A F, + gsCi A J^. (109) 

In our cases, since the forms H-^ and F4 both have overall factors of the volume-form of 
the S"^ in 123 subspace due to the isometry S0{3) x S0{6), the terms F4 AF4 and F4AH3 
are zero. The equations of motion for the fluxes on the background of the near horizon 
geometry of a shell of multi-center DO branes with D2 sources turned on are then modified 
to: 

dH^ = 0, (110) 

dGe = J7, (111) 

d[Z^\H3-gs*9G,)] = 0, (112) 
dx''Ad[-g:'Z^U,H, + {Z^'-l)G,] = gsC.AJj. (113) 
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Since J7 is a delta function located at the 5*^ with radius tq, so for the right side of 
equation (jll3p . we need to consider the Ci at ri = tq, r2 = 0. The multi-center warp-factor 
Zi on the location of the delta function is infinity, so Ci = —g~^dx^ on the location of the 
delta function, where we used the gauge choice Ci = g^^{Z^^ — l)dx^. Since dGe = J7, 
the last equation ()113|) becomes d[Zi^(*gH^ — gsGe)] = 0. 

F Solving the equation with source 

In this appendix we solve the equations in in section 4. Since *9Gq^^ is closed, it can 
be written as *9Gq^^ = {r^^dihdri + r^^d2hdr2) Auj2, so the equation dG^^^ = J7 gives 

dfh - 2r^^dih + dlh + hr:^^d2h = 2/€V^2^?;Vo'^(n - ro)S%r2). (114) 

We can convert equation (jll4j) to a Laplacian equation with source terms by making 
a derivative with respect to ri of both sides of equation ()114|1 and define the function 
Y = r^'^dih. The resulting equation becomes 

Vg^ = Q, (115) 

where Vg is the Laplacian on the 9-d flat space: Vg = ['i"i'^di{rldi) + r^^92(r|92)] and 
Q = 2K'^T£)2g~'^prQr^'^5' {ri—rQ)5^{r2) is the source term, where 5'(ri — tq) is the derivative 
of 5{ri — To) with respect to ri. Y can be solved by integration via Green's function: 

Y{'r)= f G{'r,~r')Q{'r')d^~r\ (116) 



where Gi^r , ~r') is the green function in the 9-d flat space defined as 

v9G(r,y) = (117) 
G{r,r) = (118) 

where C is a constant. 

Now let's study the y at a point (0, 0, ri, 0, 0, 0, 0, 0, r2). It is the superposition of 
all the potentials generated by the sources at points parameterized by {r[ sin 6 cos (f, 
r^sin6'sin(y9, r^cos^^, 0,0,0,0,0,0), where the first three coordinates denote xi,X2,X3 and 
the last six denote x^, ...,xg. From equation ()ll(jj) . now Y should be: 

/C 
7^2 — o -' ^ , 2 I 2^7/2 ^^^^^2^^~VoK~^'^'(^l - ro)6%r'^)d'^~r'^r[^dr[ sin 

ATTCK^TD2g;^prl f 1 1 \x'(' \a' nio^ 

' d (ri - ro)rfri, (119) 



5rir'i \ [(ri — r^)^ + r^]^/^ [(ri + r^)^ -|- r|]^/^ 

27 



where we first integrated over cP~r'2-i d6, dip. Then we use 5' {r\ — tq) = hm^[5(r^ — ''"o + 
e) — 5{r^ — — e)] and perform the integration and then take e — > 0. We then get the 
result for Y in equation ()43|1 in section 4. 

In the region where ri ^ tq, we can expand Y in terms of powers of tq: 



Y 



5ri 



4 /105ri 315rf 



3! 



,11 



' » + 5! 



8 /-4527ri 34650r3 45045rf 



,11 



+ 



We see the leading term of (fT^ is of order r^. Y = 4^g^'_ji32.g« p 



,13 



£ ^ 105ri 
3! 



,15 



0(r, 



(120) 



3r 



H. 



Tg, where /c is a constant. Then we get h = fc^^ and the leading terms of 
^^^^ and Ci^^^ are 



H. 



{II) 



Qs *9 Cf^ = Qsk 



1-3^] dn A W2- 3 



rir2 



dr2 A W2 



G. 



{II) 



gsk^[T^~3V^l 

3 

fc^[*9r3-*93V^3] 



(121) 
(122) 



We see that the leading terms of h'^^'' and G)^'' in large r region is precisely our n = —9 
solution (jlip in section 2. Since tq oc fi, they are of the third order in /i in large r region. 



^{11) 
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